GW study of the half-metallic Heusler compounds Co2MnSi and Co2FeSi 
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Quasiparticle spectra of potentially half-metallic Co2MnSi and Co2FeSi Heusler compounds have 
been calculated within the one-shot GW approximation in an all-electron framework without ad- 
justable parameters. For Co2FeSi the many-body corrections are crucial: a pseudogap opens and 
good agreement of the magnetic moment with experiment is obtained. Otherwise, however, the 
changes with respect to the density-functional-theory starting point are moderate. For both cases 
we find that photoemission and x-ray absorption spectra are well described by the calculations. 
By comparison with the GW density of states, we conclude that the Kohn-Sham eigenvalue spec- 
trum provides a reasonable approximation for the quasiparticle spectrum of the Heusler compounds 
considered in this work. 
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I. INTRODUCTION 

Heusler compounds^ attract ever-growing experimen- 
tal and theoretical attention, largely because a vast num- 
ber of such compounds have been predicted to be half- 
metallic ferromagnets (HMF), i.e., the compounds be- 
have like a metal for one spin channel and like a semi- 
conductor for the other i^"— Peculiar electronic transport 
properties are expected from such materials, e.g., huge 
magnetoresistive effects in giant and tunnel magnetore- 
sistive (GMR, TMR) devices. 

Heusler compounds are ternary intermetallic com- 
pounds with the general chemical formula X2YZ, where 
X and Y are transition-metal atoms and Z is a main 
group element. They form the cubic L2i structure (space 
group Fm3m) with a four atom basis. The half-metals 
among the Heusler compounds follow the Slater-Pauling 
rule, which connects the magnetic moment per formula 
unit m and the number of valence electrons Ny viaS. 



m — N\i 



24. 
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Most theoretical studies of these materials have been 
based on density functional theory^'- (DFT) in the Kohn- 
Sham formalism so far,^ which gives access to ground- 
state properties, such as the total energy, atomic forces, 
magnetic moments etc. It relies on a mapping of the 
real system onto a fictitious system of noninteracting 
electrons moving in an effective potential. The half- 
metallic nature found experimentally for some Heusler 
compounds is predicted correctly by DFT, together with 
a quantitative explanation of the Slater-Pauling behav- 
ior. However, it is questionable whether the Kohn-Sham 
eigenvalue spectrum can be taken as the excitation spec- 
trum of the real system. Strictly speaking, there is no 
theoretical justification for such an interpretation. In 
fact, while the band structure often resembles the exper- 
imentally determined dispersions remarkably well, there 
are important quantitative discrepancies. For example, 
the fundamental band gaps of semiconductors and in- 
sulators are usually underestimated by a factor of 2 or 



more. This raises the question if the half-metal band 
gap is also subject to this underestimation. Studies on 
Co2MnSi indicate that this is not so: the experimental 
gap is not larger than about 1 eV as inferred from tun- 
nel spectroscopy and x-ray absorption experiments 1^""— 
This value is very close to the calculated Kohn-Sham gap. 
The bandwidth of metals and the exchange splitting of 
ferromagnets are two other important spectral quantities 
which are often unsatisfactorily described by Kohn-Sham 
DFT^i^ 

There are several approaches that allow to go beyond 
Kohn-Sham DFT in this respect. For example, DFT-I-C7 
and DFT+DMFT (dynamical mean field theory in a cor- 
related subspace) employ an effective, partially screened 
interaction parameter, the Hubbard U parameter, that 
acts between electrons in the subspace of localized states 
while the rest is treated on the level of DFT.^^ The U 
parameter itself is taken in its static limit. Dynamical 
screening effects of the itinerant electrons are thus ne- 
glected. Furthermore, the Hubbard U parameter is usu- 
ally taken as an empirical parameter that is fitted to ex- 
periment, and the artificial separation into localized and 
itinerant electrons requires a double-counting correction, 
which is not uniquely defined. LDA-I-DMFT calculations 
on half-metals suggest the presence of nonquasiparticle 
states inside the half-metal gap, which may destroy the 
half-metallic character of a materialji^""— 

Another method that allows to obtain physical elec- 
tron addition and removal energies is the GW approxi- 
mation for the electronic self-energy within many-body 
perturbation theory j^^i^^ In contrast to DFT, the GW 
method is designed for spectral properties, such as the 
band structure. Typically, it opens the gap of semicon- 
ductors and insulators and gives good agreement with 
experiments. We apply this method to Co2MnSi and 
Co2FeSi, two prototypical and potentially half-metallic 
Heusler compounds, to study the effect of many-body 
corrections on their band structures. In particular, we 
will answer the question if the GW approximation in- 
creases the half-metal band gap as in the case of semi- 
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conductors and insulators or not. As already mentioned 
above, an increase may worsen the good agreement to 
experiment achieved by Kohn-Sham DFT. 

Co2MnSi and Co2FeSi are particularly interesting be- 
cause of their large magnetic moments and Curie tem- 
peratures. They are known to form the L2i struc- 
ture with a low degree of chemical disorder this al- 
lows accurate comparison between experiment and the- 
ory. The half-metallic character and integer magnetic 
moment of Co2MnSi is already predicted by DFT.™ For 
Co2FeSi, DFT calculations predict a significantly re- 
duced magnetic moment with respect to experiment and 
the Slater-Pauling value^^S DFT+U and DFT-f-DMFT 
find a magnetic moment in accordance with the Slater- 
Pauling rule and experiment with U parameters of 1.8 
and 3eV, respectively.^^'^i However, DFT+U deterio- 
rates the spectral properties of Co2FeSi compared to con- 
ventional DFT calculations. — It is the aim of this work to 
investigate to which extent many-body corrections within 
the GW method modify or confirm the predictions made 
by DFT calculations and, in particular, whether the GW 
approximation is able to rectify the magnetic moment of 
Co2FeSi without deteriorating the spectral properties. 

II. METHOD 

In this work we present one-shot GW calculations, 
which yield the quasiparticle energies E'^y. as corrections 
on the Kohn-Sham energies ej^^., 

Kk = <k + (Ck |sjc(Kk) - <c I Ck) , (2) 

where ^J^j^ are the Kohn-Sham wavefunctions and n, k 
and a are the band index, Bloch vector, and electron spin, 
respectively. The quasiparticle correction contains the 
exchange-correlation potential u^c for which we employ 
the PBE functional,— and the GW self-energy operator, 
which is given in formal notation by E^c = iG'^Wj^^ 
where G"^ and W are the Kohn-Sham Green function 
and screened Coulomb potential, respectively. The lat- 
ter is approximated by the random-phase approxima- 
tion W = v{l — vP)~^ with the polarization function 
P — —iy^, G'^G" and the bare Coulomb interaction v. 
Notably, W does not depend on spin: quasiparticles of 
both spin directions interact via the same screened po- 
tential. 

We use the fleur^'* and spex^ programs for the 
DFT and GW calculations, respectively. These codes 
are based on the highly precise all-electron full-potential 
linearized augmented-plane-wave (FLAPW) method. 
Transition-metal 3s, 3p and Si 2s, 2p semicore states are 
treated with local orbitals, although their effect on the 
spectra is small. The muffin-tin radii are set to 2.25 and 
2.31 bohr for the transition-metal atoms and Si, respec- 
tively. We employ plane-wave and angular momentum 
cutoff parameters of fcmax = 4.0bohr~^ and Z,„ax — 8. 
The DFT calculations are performed on 256 k points in 
the irreducible wedge to obtain a reliable starting point. 



The GW calculations are performed with alOxlOxlO 
k-point mesh that contains 47 points in the irreducible 
wedge with cutoff parameters for the mixed product basis 
imax — 4 and G^^^ = 3.5bohr~^, and an additional cut- 
off A/47r/wmin = 4.5 bohr""^ for the correlation part of the 
self-energy; see Ref.H^for details. We find that 50 empty 
bands are sufficient to converge the quasiparticle spectra 
to better than 0.05 eV. This is also the estimated accu- 
racy of the k-point sampling. The self-energy is evalu- 
ated with a contour integration in the complex frequency 
plane, and Eq. [2l which is nonlinear in energy, is solved 
on an energy mesh with spline interpolation between the 
points. 

The densities of states (DOS) curves are obtained with 
tetrahedron integration and convoluted with a Gaussian 
of 0.1 eV full- width at half-maximum. Binding energies 
are always taken relative to the corresponding Fermi en- 
ergy, which is determined by the condition that the DOS 
integrates to the total number of electrons from — oo to 
the Fermi energy. All calculations are based on the exper- 
imental lattice constant of 5.64 A for both compoundsi^ 



III. RESULTS 

In Fig. [U we present the PBE and GW DOS of 
Co2MnSi and Co2FeSi. In both cases, the main effect 
of the quasiparticle corrections are downshifts of the Si s 
states (between —9 and —12 eV) by 0.9 eV and the hybrid 
•p-d states (between —4 and — 8eV) by 0.8eV-0.5eV — 
see, e.g., Ref. 0|for partial DOS plots. Additionally, the 
exchange splitting of these states is reduced. 

The binding energies of the occupied d states of 
Co2MnSi remain largely unchanged. While the abso- 
lute values of the 3 d quasiparticle energies do change due 
to the exactly cancelled self-interaction error, the Fermi 
energy changes likewise so that the difference remains 
more or less the same. Close to the Fermi energy we 
find a small increase of the exchange splitting by 0.2 eV 
in Co2MnSi, which places the GW Fermi energy closer 
to the minority valence band minimum. In addition, the 
minority gap (given by the F ^ X transition) is slightly 
enhanced from 0.82 eV to 0.95 eV, and the unoccupied 
minority d states are rigidly pushed up in energy. The 
only small quasiparticle correction of the minority gap is 
noteworthy in view of the fact that semiconductor and 
insulator gaps usually increase considerably (and rightly 
so) when treated within the GW approximation. Thus, 
the apprehension that GW might worsen the agreement 
with experiment is proved wrong with this result. This 
aspect will be analyzed in more detail in the next section. 

A similar effect is encountered for Co2FeSi, but the in- 
crease of the exchange splitting and the shift of the unoc- 
cupied d states are larger than in Co2MnSi. This places 
the Fermi energy in the middle of a minority pseudogap, 
which accommodates a light band of Fe tig character. 

The quasiparticle shifts are displayed in Fig. [5J We 
see that the d states are pushed up in energy; the oc- 
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TABLE I. Magnetic moments (in /j-b), minority F — F, minority T — X, and minority spin flip gaps (in eV) of Co2MnSi and 
Co2FeSi obtained from Kohn-Sliam DFT, the GW approximation, and experiment where available. 
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FIG. 1. Kohn-Sham and GW DOS of CoaMnSi and CoaFeSi. 




-3 3 -12 -9 -6 
Kohn-Sham energy (eV) 

FIG. 2. Quasiparticle shifts as function of the Kohn-Sham 
energy. 



cupied states move closer to E-p and the unoccupied 
states away from it. Also the increase of the exchange 
splitting around the Fermi energy becomes visible. For 
Co2FeSi, the states close to the Kohn-Sham Fermi energy 
are pushed up in energy by as much as 0.85 eV. These are 



mostly of Fe d character with 25 - 50 % admixture of Co 
d character. 

Table |T] compares the magnetic moments, the minority 
r r and r ^ X transition energies, and the minority 
spin flip gaps from the Kohn-Sham and quasiparticle cal- 
culations and from experiments. The magnetic moment 
of Co2MnSi is the same in PBE and GW and matches 
the experimental value very wellj^ With the Fermi en- 
ergy located in the pseudogap, the magnetic moment of 
Co2FeSi is increased from 5.52 /zb / f.u. to 5.89 iib / f-u., 
improving the agreement with the experimental value of 
5.97^B/f-u. considerably.^^ Hence, the one-shot GW 
approach manages to correct the magnetic moment. We 
note that the orbital magnetic moment^ii^ is not taken 
into account in our calculations. 

The minority spin flip gap, i.e., the energy required 
to promote an electron from the minority valence band 
maximum to a majority state at the Fermi energy, is 
nonzero for Co2MnSi but zero for Co2FeSi due to the 
minority pseudogap. From tunnel spectroscopy of mag- 
netic tunnel junctions one deduces a spin flip gap for 
Co2MnSi between 0.25 and 0.35 Both theoreti- 

cal values are in fair agreement with these experimen- 
tal numbers. Co2FeSi does not have a spin flip gap in 
the experiment which agrees with both calculations. 
The minority F F transition energy is increased for 
Co2MnSi by 0.13 eV, whereas it essentially remains the 
same in the case of Co2FeSi. This is very different from 
the DFT+U {U = 1.8 eV) result, where the F -> F gap 
of Co2FeSi increases to 1.8 eV.— 

We compare our calculated quasiparticle spectra with 
experimental high energy x-ray photoemission spectra 
(HXPS) taken at 7.935 keV^ The fuU valence band spec- 
tra are given in Fig. [31 with the features discussed in the 
following marked by arrows. We compare only peak po- 
sitions, as a detailed analysis of the peak heights would 
require the calculation of the transition matrix elements, 
which is beyond the scope of this paper. For both mate- 
rials, the main features of the spectra are reproduced by 
the calculations. The valence band minima of Co2MnSi 
and Co2FeSi at — 12.4eV and — 12.8eV, respectively, are 
accurately reproduced by the GW calculations. Also 
the maxima of the emission from the Si s states are 
about correct. The emission maxima of the p-d hybrid 
states are in good agreement with the PBE calculation, 
whereas GW places them too low in energy compared 
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to experiment, while their onset is described better. It 
is difHcuh to assign the individual structures between 
— 5eV and E-p in the experimental spectra to the var- 
ious peaks in the quasiparticle spectra. However, the 
overall agreement seems to be reasonable in both cases. 
The plasmon frequency calculated within the random- 
phase approximation amounts to 4.7 eV and 6.0 eV for 
Co2MnSi and Co2FeSi, respectively, in agreement with 
previous calculations i^i^ These energies are well within 
the valence band region, indicating that the measured x- 
ray photoemission spectra might be affected by plasmon 
satellites. 

Additional high-resolution HXPS spectra taken close 
to the Fermi energy are shown in Fig. |4l Both spectra are 
well described by the GW calculation. For Co2MnSi, the 
main feature at —1.25 eV, arising from a Co-Mn majority 
d state, is placed 0.1 eV too high in GW. The shoulder 
at — 0.7eV arises from a pure Co minority d state and 
is reproduced by GW. The structure at — 0.3eV in the 
experimental Co2MnSi spectrum might be related to the 
minority valence band maximum, which appears at about 
the same energy in PBE and GW, see the spin-flip gap 
values in Table HI Strangely, the GW DOS does not show 
a structure at this energy in contrast to the PBE DOS. 
A comparison with Fig. [T] reveals that while the minor- 
ity DOS drops at — 0.3eV, the majority DOS happens to 
increase at exactly the same energy so as to compensate 
the decrease from the minority states. However, we note 
that even a small difference in the transition matrix ele- 
ments of spin-up and spin-down states, which have been 
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FIG. 3. Comparison of experimental high-energy x-ray pho- 
toemission spectra and total DOS of Co2MnSi and Co2FeSi. 
Experimental data taken from Ref. i27i . 
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FIG. 4. Comparison of experimental high-energy x-ray pho- 
toemission spectra and total DOS of Co2MnSi and Co2FeSi 
close to the Fermi energy. The horizontal dashed line denotes 
the additional background added to the theoretical spectra. 
Experimental data taken from Ref. [27l . 



neglected in the present work, are expected to produce a 
structure in the GW spectrum at the correct energy. 

The photoemission spectrum of Co2FeSi close to the 
Fermi energy in Fig. 2] is well described by the GW calcu- 
lation and improves on the PBE result. The features are 
less pronounced than for Co2MnSi; however, the shoulder 
at — 1.3eV and the shape of the spectrum below — 0.6eV 
are reproduced. 

Now we turn to the unoccupied states. We focus on 
the transition-metal d states, which can be mapped out 
element-specifically by soft x-ray absorption spectroscopy 
at the L3 edges (using the 2p — 3d transitions). In 
Fig. [S] we compare the experimental L3 absorption spec- 
tra of Co, Mn, and Fe in Co2MnSi and Co2FeSi with 
the corresponding GW d electron DOS. The absorption 
maxima are aligned with the DOS maxima. The shapes 
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FIG. 5. Top row: experimental spin-averaged x-ray absorp- 
tion spectra at the Co, Fe, and Mn L3 absorption edges. Bot- 
tom row: site-resolved GW d electron DOS. The absorption 
maxima are aligned with the theoretical DOS maxima. Ex- 
perimental data taken from Refs. [sO and [2^. 
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of the spectra agree with the computed DOS; also, the 
aUgnments with the Fermi energy seem reasonable, and 
the hybridizations are visible in spite of the large life- 
time broadening of the spectra. For a detailed compar- 
ison of the energy levels one would have to take into 
account the interaction of the core hole with the pho- 
toelectron, i.e., an exciton. This effect is of the order 
of 0.3. . .0.5eV, and it affects the final states in depen- 
dence on their symmetry and localization] ^^i'^^ A consis- 
tent treatment of the optical absorption process would re- 
quire solving the Bethe-Salpeter equation. Kallmayer et 
al. have taken the exciton binding energy as 0.5 eV and 
assumed the exciton to effect a rigid shift of the un- 
occupied d states towards the Fermi level. With these 
assumptions, they find that the maximum DOS of Co 
should be at 0.9 eV and 0.6 eV above E-p for Co2MnSi and 
Co2FeSi, respectively.- - These values agree with our cal- 
culated GW values within 0.1 eV, while the Kohn-Sham 
spectrum shows a larger discrepancy, see Fig. [TJ 

The unoccupied minority d states of Co2FeSi are 
mostly shifted rigidly upwards in the GW calculation. It 
was recently shown, that x-ray magnetic linear dichroism 
spectra of Co2FeSi can be described by a DFT calcula- 
tion with the PBE functional plus a rigid shift of the d 
statesi^ We conclude that the spectrum of unoccupied 
states is described correctly within the GW approxima- 
tion. 



IV. ROLE OF THE SCREENING 

In the GW approximation, the screened Coulomb in- 
teraction VF(r,r';e) is the key ingredient. Intuitively, 
one may expect that the similarity of the PBE and GW 
results arises from the metallic screening of the major- 
ity spin channel. To test this conjecture, we have com- 
puted the GW gap of Co2MnSi without metallic screen- 
ing. We also analyze the importance of local-field effects 
and briefiy discuss results from a one-shot PBEO hybrid 
functional scheme i^i^ 

Neglecting screening altogether, i.e., replacing W by 
the bare Coulomb interaction v, we obtain the (non- 
self consistent) Hartree-Fock gap of 9.65 eV, a gross over- 
estimation. Now we allow for screening effects but sup- 
press the metallic screening. We achieve this by replacing 
polarization contributions from the majority spin chan- 
nel, where metallic screening takes place, by the polar- 
ization arising from the minority spin electrons, i.e, we 
use P = 2P^. This enforces a long-range W also in the 
static limit, since the electrons cannot fiow freely in the 
gapped minority channel, which would enable them to 
screen test charges completely. Employing this artifi- 
cial semiconductor-like polarization, which exhibits a fi- 
nite dielectric constant of Eoo ~ 14, we obtain only a 
slightly larger minority energy gap of 0.97 eV. On the 
other hand, setting P = 2P^ reduces the gap to 0.86 eV. 
Clearly, the majority electrons generate a more effective 
screening, but the differences in the gap values are rel- 



atively small. Long-range metallic screening does not 
seem to contribute significantly to the total screening, 
and screening taking place at short distances seems to be 
more effective. 

To investigate this further, we exclude local-field ef- 
fects. Local- field effects arise from density fiuctuations 
of a different wave length than their generating fields. 
These couplings are related to the offdiagonal elements 
of the polarization matrix P represented in a plane- 
wave basis. (We employ, instead, a basis of eigenvectors 
of the Coulomb matrix represented in the mixed prod- 
uct basis, which are, however, reasonably close to plane 
waves.) Setting these offdiagonal elements to zero implies 
that the screened interaction W{y, r'; e) only depends on 
the difference |r — r'| rather than on the absolute po- 
sitions r and r'. This is equivalent to saying that the 
charge density within the unit cell and its screening are 
homogeneous .^93 The resulting energy gap of 1.65 eV is 
nearly twice as large as the Kohn-Sham value. Also, the 
low-lying s and p-d states are affected significantly: they 
shift by about 0.5 eV upwards in energy with respect to 
the PBE result, at odds with experiment. Furthermore, 
the exchange splitting of the occupied d states increases 
and the minority spin-fiip gap vanishes. Thus, the charge 
inhomogeneity plays a crucial role for the screening prop- 
erties. We note, that Damewood and Fong found simi- 
larly small changes of the half- metallic gaps of zincblende 
CrAs, MnAs, and MnC in the GW approximation with 
respect to PBE calculations, and a similar behaviour of 
the gap size with respect to the local-field effectsr— 

In recent years, potentials derived from hybrid func- 
tionals, e.g., PBEO,-- have often been used as an ap- 
proximation to the electronic self-energy. Being nonlocal 
they fulfill an important condition of the self-energy. Hy- 
brid functionals have been shown to overcome the typical 
underestimation of band gaps within Kohn-Sham DFT. 
However, dynamical effects are not taken into account, 
and screening is only considered in an average way by 
the parameter that mixes the nonlocal and local parts. 
In the PBEO functional this mixing parameter is univer- 
sally taken to be 0.25.^^^ Since the GW approximation 
contains the bare exchange exactly, we can easily calcu- 
late a one-shot (non-selfconsistent) PBEO energy spec- 
trum. We find that while PBEO gives similar results for 
the binding energies of the low-lying s and p-d states as 
the GW approximation, it completely fails in determin- 
ing the minority gap, for which it yields 3.03 eV. Also, 
the exchange splitting is strongly overestimated in both 
cases. Thus, only a dynamical self-energy can simulta- 
neously describe states close to the Fermi energy and far 
away equally well. 



V. CONCLUSIONS 

We have presented one-shot GW calculations of the 
(potentially) half-metallic Heusler compounds Co2MnSi 
and Co2FeSi. The GW quasiparticle spectra are qualita- 
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tively similar to the Kohn-Sham eigenvalue spectra, but 
show important quantitative differences. In particular, 
the GW approximation predicts an electronic structure 
with a minority pseudogap in the case of Co2FeSi, which 
corrects the magnetic moment per unit cell to nearly an 
integral number, consistent with available experimental 
data. 

The quasiparticle spectra are in good agreement with 
photoemission and x-ray absorption data for both com- 
pounds. The electronic screening is effective at short dis- 
tances and charge inhomogeneities play an important role 
for the screening. Furthermore, it has been shown that 
the PBEO hybrid potential cannot be used as an approx- 
imate self-energy: it even yields worse results than the 
local PBE potential. 



So far, most theoretical studies of Heusler compounds 
have been based on the Kohn-Sham band structure. In 
this work, we have demonstrated that it can, in fact, 
represent a reasonable approximation to the many-body 
quasiparticle spectrum, which confirms previous success- 
ful calculations of spectral properties of Heusler com- 
pounds within DFT. 
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